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We propose and analyze a scheme to observe topological phenomena with ions in microtraps. We 
consider a set of trapped ions forming a regular structure in two spatial dimensions and interacting 
with lasers. We find phonon bands with non-trivial topological properties, which are caused by the 
breaking of time reversal symmetry induced by the lasers. We investigate the appearance of edge 
modes, as well as their robustness against perturbations. Long-range hopping of phonons caused 
by the Coulomb interaction gives rise to flat bands which, together with induced phonon-phonon 
interactions, can be used to produce and explore strongly correlated states. Furthermore, some of 
these ideas can also be implemented with cold atoms in optical lattices. 
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I. INTRODUCTION 

Atomic systems have revealed themselves as ideal labo- 
ratories to study many-body quantum systems. Atoms or 
molecules in optical lattices [l|, Q realize Hubbard mod- 
els, both bosonic [3] and fermionic and thus can be 
used to study some of the most basic models of electron 
transport and magnetism appearing in condensed mat- 
ter physics. Starting from the seminal experiment [Bl, 
many groups are now able to implement different ver- 
sions of Hubbard-like models. Although at the moment 
the temperatures in most of those experiments are too 
high, single site addressing [6, 7] and other techniques 
are currently opening up new avenues for the observation 
of many new and exciting phenomena. Trapped ions can 
also realize bosonic Hubbard models \§\ , as well as many 
other spin models [ol, |lo|. Despite recent experimental 
breakthroughs [ll| , ions are not yet as versatile as atoms 
in lattices. However, this situation may soon change with 
the development of surface traps |12|-[l6j or through the 
combination of ions and optical lattices [l^ • 

One of the goals of atomic experiments is the simula- 
tion of the Fractional Quantum Hall Effect (FQHE), one 
of the most celebrated discoveries in condensed matter 
physics. In the standard electronic scenario, it appears at 
very low temperatures in the presence of both an external 
magnetic field and strong interactions. The first restricts 
the electron motion to the first Landau level, whereas the 
latter is responsible for the generation of strong correla- 
tions. It has recently been suggested that, under some 
conditions, in lattice systems a similar behavior should 
occur even in the absence of a net magnetic flux [l8|, [l9| . 
The first condition is that the band structure should pos- 
sess a non-trivial topology, as characterized by the so- 
called Chern Number (CN). The second one is that the 
width of the band should be much smaller than the band 
gap, i.e., one has flat bands, something which is asso- 
ciated to hopping beyond nearest neighboring sites. In 
such a situation, the lattice problem resembles the stan- 
dard scenario for the appearance of the FQHE, whereby 
the flat band plays the role of the first Landau level. In- 



teraction energies larger than single particle ones should 
thus give rise to strongly correlated states at low temper- 
atures and, in particular, to states displaying fractional- 
ization as well as to anyonic excitations. 

Even in the absence of interactions, the non-trivial 
topology of the energy bands gives rise to very intriguing 
behavior. In particular, when we consider a problem with 
open boundary conditions (as opposed to periodic bound- 
ary conditions) new energies appear within the band gap 
in the spectrum, with a nearly linear dispersion relation. 
The corresponding modes are located at the boundary of 
the system, and thus those are called "edge modes" [20| . 
The presence of those modes is presumably very robust 
against perturbations [2l|, [22| . 

In atomic systems, in the bulk, one can simulate the 
action of a magnetic field either by rotating the trap hold- 
ing the atoms [23] or by using laser fields [2J-[26|. In 
optical lattices, a very promising alternative consists of 
using complex hopping amplitudes induced by the lasers 
[2?! (see also [26|). The appearance of FQHE-like behav- 
ior in optical lattices was first considered by Hafezi et al 
[28| in the low density limit, where the lattice does not 
play an important role. In addition, they proposed to 
use dipole-dipole interactions to increase the band gap 
(see also [22, 29-31]), and thus be able to increase the 
density. In a recent paper, an implementation based on 
polar molecules [32] is proposed which achieves very flat 
bands by exploiting the dipole-dipole interaction among 
them, and thus is very well suited for the observation of 
the FQHE. 

In this paper we show that trapped ions can be used 
to investigate some of the intriguing features mentioned 
above. We consider ions in microtraps, forming a regu- 
lar structure in two spatial dimensions. They are driven 
by lasers, which are far off resonance with respect to 
an internal electronic transition. By carefully choos- 
ing the laser configuration, the ions experience a force 
that breaks the time-reversal symmetry. Furthermore, 
phonons (i.e., motional quanta in each of the ions) can 
hop from one ion to another due to the Coulomb interac- 
tion, which has a long range. As a result of all that. 
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we find a band structure with a non-trivial topologi- 
cal structured, as witnessed by the corresponding CN. 
We investigate the appearance of edge modes, as well 
as their fragility in the presence of perturbations. We 
also study the regime of parameters where the bands are 
very flat, whereby FQHE-like behavior could be investi- 
gated by adding phonon-phonon interactions as proposed 
in [8.]. Note that other methods to obtain flat bands (al- 
beit without topological properties) have been recently 
proposed for trapped ions using time-dependent gadgets 
[335. We emphasize that the scheme presented here is 
quite different from the one proposed for atoms in optical 
lattices, as here the hopping terms are all real, whereas 
the on-site Hamiltonian is the one that breaks the time- 
reversal symmetry. In fact, our scheme resembles more 
that studied in the context of photonic edge modes [34^. 
The configuration we propose here could also be suitable 
in setups involving atoms in optical lattices. 

This paper is organized as follows: In Sec. II, we intro- 
duce the model setup and derive an effective Hamiltonian 
for the phonons. In Sec. Ill, we analyze the energy spec- 
tra of the system with toroidal, cylindrical, and planar 
geometries. In particular, we derive the CN for the dif- 
ferent bands and relate them to the appearance of edge 
modes. In Sec. IV, we propose a scheme to detect the 
edge modes by using a periodically driven electric field 
that shakes the ion lattice. In Sec. V, we discuss that, 
due to the flat lowest band, phonon-phonon interactions 
could lead to the observation of a bosonic version of the 
FQHE. Some of the details of the derivations are pre- 
sented in the appendices. 



II. MODEL SETUP 

We consider N ions confined in microtraps and in- 
teracting with lasers in a planar geometry. The lasers 
are tuned far off-resonance with respect to any optical 
transition, so that the ion stays most of the time in its 
ground electronic internal state. Virtual transitions me- 
diated by the lasers modify the atomic motion. In par- 
ticular, they give rise to a Doppler effect which, when 
combined with dipole forces, generates something that 
resembles the Lorenz force produced by a magnetic field. 
This force, which breaks the time-reversal symmetry of 
the problem, is one of the main ingredients required for 
the appearance of the topological features which we are 
looking after. The other one is the Coulomb interaction 
among the ions, which allows phonons to move from ion 
to ion [8]. The effective dynamics of the phonons in the 
ion lattice are described by a Hamiltonian which breaks 
time-reversal symmetry and, as we will see, has the de- 
sired topological properties. 

There are many lattice geometries that give rise to 
topological features. Here, for the sake of concreteness, 
we will consider a hexagonal (honeycomb) geometry [see 
Figdl^a)]. In this section we will first introduce the in- 
teraction of the laser with each ion, and derive the cor- 
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FIG. 1: (Color online) Configuration: (a) Ions are trapped in 
microtraps forming a planar hexagonal geometry, (b) Laser 
and magnetic field configuration for each ion in a micro- 
trap. There is a homogenous magnetic field as well as a laser 
standing-wave with circular polarization along the x direc- 
tion, and two counter-propagating lasers in lin±lin configura- 
tion along the y direction, (c) Ion internal structure, where 
the quantization axis is chosen along the magnetic field (x 
axis). For simplicity, we consider no nuclear spin (hyperfine 
structure). The magnetic field generates a Zeeman shift. The 
lasers propagating along the x direction (in red) produce dif- 
ferent AC- Stark shifts on the magnetic levels. Lasers along 
the y direction induce Raman transitions, (d) Detail of the 
hexagonal lattice, with the unit vectors ex,y indicating the two 
oscillation directions of the ions (vibrational modes). Differ- 
ent colors are used for the A and B sublattices. 



responding Lorenz-like force. Some of the details regard- 
ing a particular laser configuration that implements those 
forces are given in Appendix A. Later on, we will con- 
sider the interaction among the ions in order to derive the 
effective Hamiltonian describing the phonon dynamics. 



A. Single-ion dynamics 

Let us first consider a single ion, as shown in Fig. [T]3, 
with two electronic (internal) ground levels 1^1,2) inter- 
acting with lasers. The single-ion Hamiltonian contains 
three parts: H^i = Huo + H^^^ + H^y^^ where {H = 1) 
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Here, i^no describes the harmonic oscihation of the ion 
in the plane, where M is the mass, and Tq,, Po^^ and ujo^ 
denote the position and momentum operators, and the 
trap frequency along the direction a, respectively. The 
terms H]^'^ describe internal level-dependent dipole and 
radiation pressure forces, respectively. Here K is twice 
the wavevector of the lasers, ^x,y the effective Rabi fre- 
quencies, cr^ = 1^2) {g2 \ - \9i) and cr+ = \g2) {gi\. 

Hamiltonian ([1]) can be implemented in different ways. 
One possible one is depicted in Fig. mb,c), which uses a 
magnetic field along the x direction, and two pairs of far 
off-resonant counterpropagating beams along the x and 
y direction. In Appendix A we give the details of this 
scheme, as well as how to obtain H^i after the adiabatic 
elimination of the excited states. 

It is convenient to transform the Hamiltonian ([1]) ac- 
cording to the unitary operator U = e-x.-p{iKry(Jz /2) 



KPy, 



U^H^iU = Hno + ^cT, + {n.Kr,+ '-^)a,, (2) 



where the last term represents the Doppler shift along the 
y direction and we have used the standard notation ax = 
(7+ + <j~ . In the following we will simplify the model by 
assuming that \^y i ^ Vx, ^yVy^ where r]a = Ki^ 
are the dimensionless Lamb-Dicke parameters, which are 
typically < 1 (and ia = (2Mco'q;)~^/^) is the ground state 
size). Assuming that the ion is initially prepared in the 
state \D = (1^2) — |^i))/>/2, we show in Appendix B that, 
after the adiabatic elimination of the internal states, one 
ends up with the effective Hamiltonian 

i^eff = ^HO + ^rxPy. (3) 
The position-momentum coupling is characterized by 



■E 

a=x,y 



(4) 



and ^HO describes harmonic oscillations with renormal- 
ized frequencies uja (see Appendix B). This Hamiltonian 
resembles the one describing the Lorenz force exerted by 
a magnetic field on a charge (note, however, that it is not 
invariant under rotations along the z axis). 



B. Coulomb Interactions 

Let us now consider the whole set of ions. As men- 
tioned above, the microtraps are disposed according to a 
2D honeycomb lattice, as shown in Fig. [Ha,d). The in- 
teraction of each of the atoms with the lasers is given by 
(j3|). Coulomb forces induce repulsion between the ions. 
Assuming that d £a^ we can expand the Coulomb in- 
teraction up to second order in the displacements ri*^ 
about the center of the i-th microtrap, obtaining 



a,l3' p 



(5) 



where 

(6) 

As shown in Fig. [l]i, Gq, are the unit vectors, and Yiijd 
is the vector pointing from the site i to the site j. We 
have omitted the linear terms in ra\ which just redefine 
the trap center. Furthermore, the /7^'*^ just shift the 
frequencies cD^,, and can be absorbed in cDq,. 

The complete Hamiltonian, consisting of the single- 
ion terms as well as the Coulomb interaction, can be 
rewritten in terms of the phonon annihilation operators 
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2McD, 



(7) 



In the stiff limit j3cx = / {MCbl^d^) <C 1, where we can 
drop the non-phonon number conserving terms [9], and 
we finally obtain the Hamiltonian for the phonons 



(8) 



This describes phonon hopping without interaction (i.e., 
it is quadratic in creation and annihilation operators), 
and is characterized by the "single-particle" Hamiltonian 



a, (3 



27a 7/5 



• ^ij{^ax^f3x + Jy^ay^fdy) 



-iVhSij{5ax^(3y - ^ay^(3x) 



(9) 



in terms of the dimensionless parameters Px-, — 
\/Coa/Cbx^ and Vb = ^72y/(2a)x)- For the sake of sim- 
pHcity, in the following we will set 7^ = 1, so that the 
problem is fully determined by two parameters. The cor- 
responding values are restricted by the conditions em- 
ployed in the derivation of the Hamiltonian ([9]). In 
particular, /^^^ <C 1, and l^^l = rj'^Qx/^y ^ ^7 where 
we have taken rjx = r]y = rj. Taking the ^^Ca+ion, 
uox = ujy = Vtx ^ O.lMHz, Qy ^ 0.3-lMHz, and 
d ~ 10-100/im, we obtain values for 77 up to 0.8 |35|, 
Px e [0.01,0.3] and \Vb\ G [0.05,0.3], which are the val- 
ues that we take in the examples given below. 

The last term in Eq. (|9]) breaks the time-reversal sym- 
metry, which results from the position-momentum cou- 
pling. Thus, a non-trivial topological phase may be ex- 
pected to emerge in this setup. 



III. TOPOLOGICAL PROPERTIES 

In this section, we study the topological properties 
of the ion system. For that, we first consider periodic 
boundary conditions (PBC), i.e., a toroidal geometry, 
and analyze the spectrum of the Hamiltonian. Given 
the periodicity of the problem, it has a band structure. 
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FIG. 2: (Color online) Different configurations considered in 
the text: (a) torus (PBC); (b) cylinder. 

where the existence of a non-trivial topological structure 
of each band is characterized by the CN [20]. We find 
that, depending on the values of and V5, some of the 
bands exhibit non-zero CN. Then, we consider a cylindri- 
cal geometry (i.e. PBC along one direction only). The 
spectrum is very similar to the one corresponding to the 
torus, but now new energies appear at the gaps between 
bands with certain CN configuration. The correspond- 
ing eigenmodes (so-called edge modes [20|, [36|) appear in 
pairs, are localized at the two boundaries of the cylinder, 
and propagate in different directions. Finally, we con- 
sider the experimentally relevant situation of full open 
boundary conditions. In such a case, we also find edge 
modes. We further investigate their robustness with re- 
spect to small random perturbations. 

A. Band structure and Chern numbers 




FIG. 3: (Color online) Single-particle energy spectra in units 
of Cox for the system with PBC and N — 400. The left panels 
show the energy spectra as a function of the quasimomenta kx 
and ky (in the first Brillouin zone), whereas the right panels 
display the same spectra but as a function of kx only, in order 
to highlight the different bands and the corresponding gaps. 
The Chern numbers of each band are shown to the right of 
such plots, (a) Px = 0.02 and Vb = -0.1; (b) Px = 0.02 and 
Vb = -0.2; (c) Px = 0.04 and H = -0.1; (d) Px = 0.04 and 
14 = -0.2. 



In this subsection, we study the topological properties 
of the system with PBC, namely, the system has transla- 
tional invariance along both x and y directions [see Fig. 
[2](a)]. In order to find the single-particle spectrum of 
H (i.e. that of 1-L) we use the Fourier modes, and thus 
diagonalize 

^".'3(k) = l^^^>-*'^-^-, (10) 

where is the number of sites and k = {k^^ky) is 
the two-component quasimomentum in the first Brillouin 
zone (BZ). 

The energy spectrum of ^^^/^(k) is shown in Fig. [3] for 
different parameters j3x and Vb- As one can see in the 
plot, there are a number of energy bands. For /^^^ = 0.02 
and V\y = —0.1, there is a big gap between the first and 
the second bands, and a small one between the second 
and the third bands. For /3x = 0.02 and Vb = —0.2, the 
gap between the second and the third bands becomes 
larger, and the gap between the first and the second one 
remains open. For larger /3x ^ 0.04 and Vb = —0.1, the 
band structure remains the same. However, for /^^^ ~ 0.04 



and Vb = —0.2, the gap between the second and the third 
bands becomes closed, whereas that between the first and 
the second band remains. 

The topological properties of each band, 5, are char- 
acterized by the CN, which is defined in terms of the 
eigenmodes of |5, k), in the limit A/" ^ 00 as [20| 

Cs = / — Vk X A,.. (11) 
Jbz ^tt 

Here, the integration is taken over quasimomenta in the 
first BZ, and Ai^ = —i (s, k| 9k, s |<§, k). The CN can be in- 
terpreted as the Berry phase obtained if we change adia- 
batically the quasimomentum cyclically along the bound- 
ary of the BZ [37]. Alternatively, if we interpret Ai^ as a 
vector potential, the C is nothing but the total number 
of magnetic fluxes through the whole BZ. Because the 
Brillouin zone forms a closed surface, this number (and 
thus the CN) must be an integer, as it is required by the 
single- valuedness of gauge transformations [38[. 

Given that the CN is an integer, we do not need to 
calculate it for A/" ^ oc, but a finite N can already give us 
its value with very high confidence. In general, however. 
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FIG. 4: (Color online) As in Fig|3]for the cylinder geometry 
and N = 400: (a) /^^ = 0.02 and Vb = -0.1; (b) /3x = 0.02 
and Vb = -0.2; (c) = 0.04 and Vb = -0.1; (d) = 0.04 
and Vb = —0.2. The edge modes are marked in blue. 



it is very hard to determine the CN through Eq. ([TT]) 
numerically, since Ai^ contains singularities as long as 
C 7^ 0. Thus, we use the alternative expression [39| 



= Im / d^k 



E 

Sl<S<S'^ 



(12) 



where Si = 1,2,... denote the bands from the bottom to 
the top, (Sfe,Hk)si,s2 = (si,k| S/c,Hk |52,k), and £^k,Si is 
the energy corresponding to k). Starting from Cg, we 
can recover the CN using 



(13) 



We have determined the CN of the three bands ap- 
pearing in [3l For a given as |Vb| increases, the CN 
of the first, second and third bands vary from (—1,0, 1) 
to (—1,1,0). A non-zero CN indicates the presence of 
a non-trivial topological property in the system which, 
as we will see in the next subsection, is reflected in the 
appearance of edge modes. 



B. Edge modes 

Based on the non-vanishing CN found in the previous 
section, we expect the emergence of edge modes in the 
boundaries of our system. In this subsection, we focus 
on those modes for two different geometries: cylinder 
[2fb)] (where there are PBC along the x-direction only), 
and plane [Ha)]. 

For a cylindrical geometry, the translational invariance 
along the y direction is broken. Still, given the periodicity 
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FIG. 5: (Color online) As in Fig |3] for the planar geometry 
and 384 sites: (a) /^^ = 0.02 and Vb = -0.2; (b) /3a: = 0.04 
and Vb = —0.1. Here, / denotes the eigenmode. Note that for 
384 sites, there are 768 eigenmodes. 



along X, we can use quasimomenta in this direction, and 
thus consider the single-particle Hamiltonian 



1 



ikx (^ix jx ) 



(14) 



where 



denote the x and y coordinates of ion 



kx = 27rn/Nx with n < Nx integer, and Nx are the 
number of sites in a row. By numerical diagonalizing 
(^x) foi" each value of kx we obtain the energy spec- 
tra shown in Fig. |H They are to be compared with those 
of Fig. [3l We choose two typical examples to explain the 
relation between the CN and the edge modes. We first 
consider f3x = 0.02 and Vb = —0.2. The corresponding 
spectra are shown in Fig. [3jb) and Fig. HJ^b), respec- 
tively. As explained in the previous section, the CN for 
the first, second, and third bands are —1,1, and 0, respec- 
tively. Thus, in the cylinder, two edge modes located in 
the up- and bottom- boundaries emerge in the first band 
gap. They are signaled by the blue dots that display 
a nearly Hnear dispersion relation with opposite slopes, 
which indicate that they are counter-propagating. If we 
inspect the band gap between the second and the third 
band, we see that no edge modes show up. The reason is 
that total CN of the first and the second bands (i.e., its 
sum) is zero. The second example is for ^x = 0.04 and 
Vb = -0.1 [Fig. IH^c)]. Here, the CN for the first, second, 
and third bands are —1, 0, and 1, respectively. There- 
fore, in each band gap there are two counter-propagating 
edge modes. 

For the planar geometry, the translational invariance 
is completely broken. By numerical diagonalizing the 
Hamiltonian 



E 

3,13 



(15) 



we obtain the spectra of Fig. [5] and the corresponding 
eigenmodes, uK As shown in Fig. [5fa), for /3x = 0.02 and 
Vb = —0.2, a set of edge modes emerges in the (small) 
first band gap, while in the large second band gap there 
are none. Figure [51(b) shows that edge modes emerge 
in both the first and second band gap for /3x = 0.04 
and Vb = —0.1. Again, these results can be understood 
in terms of the CN corresponding to the three bands. 
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FIG. 6: (Color online) Spatial distribution for typical edge 
and bulk eigenmodes corresponding to the system with planar 
geometry, with 384 sites, Px = 0.04 and Vh = —0.1. (a,c) 
For an energy E = O.TcDx, corresponding to an edge mode 
and (c,d) E — O.QcJx, corresponding to a bulk mode; (a,c) x 
component (b,d) y component. 



FIG. 7: (Color online) Single particle spectra in units of ujx 
(a,b) and the spatial distribution for the x and y components 
(c,d) of an eigenmode with energy E — O.Ta;^ for a planar 
geometry with 384 lattice sites. We have added to the values 
Vh and I5x random numbers as explained in the main text. 



In order to verify that those modes are located at the 
boundaries of our system, we plot in Fig. [6] the spatial 
distribution of a typical eigenmode corresponding to the 
edge (a,b) and bulk (c,d). That is, for each ion i, we plot 
at its location [i.e., at coordinates {ix-,'iy)] the value of 
for a = X (a,c) and a = y (b,d). 
Edge modes are supposed to be robust against small 
perturbations (with energy densities smaller than the 
band gap; or, strictly speaking, as long as the gap does 
not close) [2l|. One can understand this behavior in 
a simple way by resorting to the cylindrical geometry. 
There, impurities (or perturbations) cannot back-scatter 
the propagation of phonons in one boundary since there 
are no modes available (they are located in the other 
boundary of the cylinder). In order to test this robust- 
ness, we have introduced random perturbations in the 
planar geometry. In particular, we have added for each 
of the ions: (i) to l^^l a random number in the interval 
(0.15,0.25) and to Cox a random number in the interval 
(0.95, 1.05) in Fig. EJa); (ii) to IV5I a random number in 
the interval (0.05,0.15) and to uJx a random number in 
the interval (0.95, 1.05) in Fig. [7l^b)-(d), which results in 
a random variation in Px- As we can see, in Fig. [71^a,b) 
the energy spectrum still has some energies within the 
band gap, and the corresponding modes are located at 
the boundary [Fig. [7l^c,d)]. This confirms that the edge 
modes are well-protected and robust. 



IV. PROBING THE EDGE MODES 

In this section, we propose a method to detect the edge 
modes. Since the ions are charged, they can respond to 



an external electric field. We consider a periodic electric 
field that shakes the lattice along the x-direction with fre- 
quency ujd. If we initially prepare the ions in the motional 
ground state (i.e., the state with no phonons), phonons 
will be generated by the electric field. After a time tf , the 
field is turned off. As we will show below, if the frequency 
Ud is chosen to be the one corresponding to the edge 
modes in the spectrum, phonons will be generated at the 
boundary of the system, indicating the presence of edge 
modes at that frequency. The presence of phonons can 
be detected using standard ion trap techniques [40|, |4l| . 

The interaction between the ions and the electric field 
reads 



(16) 



where Fj is the force driving the ion at site j. We de- 
fine the density distribution of the x-phonon at time t as 
Pj{t) = {0\ai^^\t)ai^\t)\0). Here, |0) is the (vacuum) 

state with no phonons, and ax\t) denote the operators 
in the Heisenberg picture. By solving the Heisenberg 
equations of motion, we obtain 



E 



-iEit 



El + mcjd 



(17) 



where fi = QjU^j^^ is the effective driving force corre- 
sponding to the eigenmode Qj = Fj j \j2Mujx' 

In Fig. [HI we plot the normalized phonon density dis- 
tribution 



Pj(^final) = 



Pj(^final) 



(18) 
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FIG. 8: (Color online) Density distribution of x-phonons after 
the evolution time tf for px = 0.04 and Vb = —0.1: (a) ujd = 
2.{)Cjx] (b) cjd = 0.60)0;; (c) cjd = l.OcJa^; (d) uj^ = O.TcJa,. 



for different driving frequencies uo^ at time tf = lOOO/cDa^. 
Here, we have considered the simplest case, namely, 
Vtj = Qq^ the same for all ions. In Fig. [Hl^a), the driv- 
ing frequency uJd = '^-OoJx lies outside the whole energy 
bands, and thus no resonance occurs, which results in 
a weak phonon excitation. When the driving frequency 
UJd = 0.6uJx is on resonance with one of the bulk states, 
many phonons are generated in the bulk, as shown in Fig. 
[8](b). For driving frequencies = oJx and O.YuJx some of 
the edge modes become resonant, and thus the density 
of phonons in the boundary is much larger than that in 
the bulk, which indicates the presence of an edge mode. 



V. PHONON-PHONON INTERACTION 

As we showed in previous sections, with trapped ions 
it should be possible to obtain bands in the spectrum 
with non trivial CN. Furthermore, due to the nature of 
Coulomb interactions, phonons can hop over long dis- 
tances, which may give rise to fiat bands. Thus, by en- 
gineering interactions among the phonons one should be 
able to induce strong correlations and investigate phe- 
nomena like fractionalization with this bosonic system. 
In this section we will analyze the flatness of the bands 
for the model introduced in Sec. II, as well as how to 
induce interactions among the phonons. 

The flatness, F, of a band can be characterized by the 
ratio between the band gap and the bandwidth. In Fig. 
[91(a) we have plotted F as a function of Vb and /3x for 
the first band which has CN —1. This figure shows that 
values of the order of 10 can be reached within our setup. 
Those values should be large enough to allow for the in- 
vestigation of strong correlation effects with this system. 




i^^i 0.3 i^^i 0.3 

(a) (b) 

FIG. 9: (Color online) (a) Flatness F of the first band as a 
function of Vb and I3x. (b) Same for a model where phonon 
hopping is restricted to nearest neighbors. 

In Fig. [91(b) we plot F for the same model as the one con- 
sidered here, but where phonon hopping is (artificially) 
restricted to nearest neighbors. As one can see, the flat- 
ness of the band becomes smaller than in the previous 
case, indicating that the long range hopping present in 
the ion trap system helps to flatten the band. Note that 
the model with NN hopping could be implemented using 
atoms in optical lattices and, despite the fact that the 
flatness is smaller, it may lead to interesting regimes for 
that system. 

In Ref. [8] it was shown how a laser can induce phonon- 
phonon interactions in the context of the Hubbard model. 
One can use the same idea in our setup. In particular, in 
Appendix C we show how to implement the interaction 
Hamiltonian 

Fu = C/i„tE«i''^^'4^'^'' (19) 

j 

where Uint is determined by the Rabi frequency of the 
standing wave and the Lamb-Dicke parameter. In or- 
der to enter the interesting regime of strong correlations, 
the parameter Uint should be smaller than the band gap, 
so that the interactions do not mix bands. However, it 
should be larger than the bandwidth, so that the inter- 
actions dominate over single phonon dynamics. Taking 
Cjx = O.lMHz, Px = 0.04 and 14 = -0.2, Fig. [S^d) shows 
that Uint — lOKHz should fulfill the conditions to obtain 
strong correlations. 

VI. CONCLUSIONS AND OUTLOOK 

In this paper, we have proposed a scheme to induce 
topological effects in a system of ions trapped in micro- 
traps. The main idea is to use laser fields that break 
time-reversal symmetry. In particular, we have consid- 
ered a pair of counterpropagating lasers which are cou- 
pled to the position and momentum of the ions, respec- 
tively, giving rise to a Lorenz-like force. This coupling is 
very similar in spirit to the one proposed and used in neu- 
tral atom systems. In a honeycomb lattice, one obtains 
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a band structure with non-trivial CN, which give rise to 
edge modes. We have analyzed the robustness of such 
modes against small perturbations, an proposed a way 
for detecting them. We have also studied the flatness 
of the bands, and show how to induce phonon-phonon 
interactions. Our results indicate that as soon as micro- 
traps are available, they should be very well suited to 
investigate the physics of the FQHE with them. We also 
emphasize that the setup used here to induce the Lorenz- 
like force could be also used in the context of atoms in 
optical lattices. 



(along the x-direction) induces transitions \gi) ^ \ei) 
and 1^2) ^ 1^2), according to 



1 
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\e2){g2\)e-'^'y-'^^'^}i.c.. 

(A3) 

Note that we have used that the wavevectors of the stand- 
ing and travelling waves are approximately the same, Kq. 

In the large detuning limit, where all the detunings are 
much larger than the Rabi frequencies, by adiabatically 
eliminating the six higher energy levels, we obtain the 
effective Hamiltonian 
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sin {Kqx - 
+%|52>' 



(/.)(Ai \gi){gi\ -A2 |<72)(<72|) 

H.c, (A4) 



where K = 2Ko, = n]:;g^nf^gJ[2{L0R - L02)], and 



Appendix A: Physical Implementation 

In this appendix, we analyze a possible implementa- 
tion of our model ([T]). We consider a single earth al- 
kali ion in a microtrap, in the presence of a homoge- 
neous magnetic field along the x direction, as shown in 
Fig. [TJ For simplicity, we ignore the nuclear spin, so 
that the ion has two electronic ground state levels, \gn) 
{n = 1,2), and two sets of excited levels, \en) (n = 1, 2), 
and \en) {n = 1, ... ,4), corresponding to the fine split- 
ting. We have chosen the quantization axis along the 
direction of the magnetic field (x-direction); so, for in- 
stance, \gi,2) = \ J = 1/2, mj = ^1/2)^. We set the en- 
ergy of the ground states as and ujz (Zeeman splitting), 
and that of the excited states as cjn, and ujn- 

The ions interact with two sets of counterpropagat- 
ing beams along the x and y directions. The first ones 
form a standing wave with frequency cjgw and circular 
polarization cr+, so that the electric field is Eq sin(i^rox + 
0) cos(co'sw^). They induce transitions \gi) ^ |e2), les), 
and 1^2) ^ 1 64), which are described by the Hamiltonian 



i^sw = [^63^1 (^) les) {gi\ - 

+ ^64^2 (^) 1^4) (^2 I 



^e2gi{x) |e2) {gi\ 

f H.c] cos(a;sw^), (Al) 



where ftij{x) = ftij sin(iiro^ + 0), and the Rabi frequen- 
cies Qij = (i| d • Eo \j) are determined by the amplitude 
of the electric field and the dipole moment d of ions, and 
X is the position operator of the atom in the x direction. 

The right-moving laser plane wave (propagating along 
the y direction) is vertically polarized (along 2:), has 
frequency cjr, and induces transitions \gi) ^ |e2) and 
1^2) ^ according to 



2^ ' 



R 

^29 



,\e2){gi\ + nl^Je^){g2\)e 



iKoy — iuJnt _ 



H.c. 

(A2) 

where Qfj are the Rabi frequencies, and y is the position 
operator of the atom in the y direction. The left-moving 
plane wave with frequency cjl and horizontal polarization 



S = UJR 



■ c^z - 

■ (^L 



I 

r^e2gi I 



4(CJR - UJ2) 



(A5) 



such that we can neglect all high-oscillating terms. The 
AC stark shifts 



Al 



l^esgi I 



■ 002 OUs 



\a 



646^2 I 



4(0)4 - - ^sw) 



(A6) 



can be taken to be the same, i.e., Ai = A2 = A, if we 
properly choose the frequency uj^w Taking the Lamb- 
Dicke limit, where the ion displacements Tx and Ty along 
the X and y directions are much smaller than the optical 
wavelength, and assuming that the trap center coincides 
with the node of the standing wave {(j) = — 7r/4) and 
S = —A, the effective Hamiltonian becomes 



(a+ 



AKr 



H.c.) 



las 



las 



(A7) 



where the effective Rabi frequency Q.x = A/2. Notice 
that the lattice spacing, (i, fulfills K^d = 47rn, where n 
is an integer. This ensures that each ion experiences the 
same AC-Stark shift produced by the lasers propagat- 
ing in the x-direction (cf. H]^^). For the second term, 
i.e., the Hamiltonian H^^^ along ^/-direction, we replace y 
by the displacement r^y, since the phase e*^^° generated 
by the coordinate yo of the equilibrium position can be 
absorbed into the definition of cr+. 



Appendix B: Adiabatic elimination of the internal 
states 

In this Appendix, we show the details of the adiabatic 
elimination of the internal states leading to Eq. (|3]). Fol- 
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lowing the Frohlich transformation [42|, we use the gen- 
erator S to rewrite the Hamiltonian 



Appendix C: Induction of phonon-phonon 
interactions 



^eff = e '^^singlee'^ 

= -f^single + [-f^single? + ^ I [^single 7 5*] , 5*] + • 

= Ho + {Hi + [Ho, S]) + [Hi + ^[Ho, S],S] 



(Bl) 



where Hq = Huo + Qyaz/2 and Hj = —K^Qx'^x + 
Py/2M)(Jx' In order to make the notation more trans- 
parent, we redefine the Pauh spin operators as fohows: 
(^x and CFz —CTx- By requiring 



Hi + [Ho,S]=0, 
we obtain the generator 



(B2) 



In this Appendix, we show how one can induce phonon- 
phonon interactions by using a laser. We basically follow 
the proposal of [8] . 

An additional laser standing wave propagating along x- 
direction with frequency cjgw and a~ polarization induces 
the transitions \gi) ^ \ei) and \g2) ^ |ei), |e2). The 
light-ion interaction is 

iJint = [^eigi{x)\ei) {gi\^neig2{x)\ei) {92] 

+1^62^2 (^) 1^2) {g2 \ +H.c.)cos(a;sw^), (CI) 
where the Rabi frequencies 

flij{x) = Ctij cos{Kox)^ (C2) 



s = n^n^T^^ + TT^) 

ily — UJX ily -\- UJx 

Here, defined according to Eq. ([7]). In the limit 

\Vty^uJa\ ^ ^xVx^ ^yVy^ we can keep up to second order 
in the expansion (|Bip and obtain the effective Hamilto- 
nian 



(B4) 



In that limit, if initially the internal state is the eigenstate 
of cr^, we can project the Hamiltonian i^eff onto that 
eigenstate and obtain 



pI 



1 



xPy-) 

(B5) 



where the renormalized frequencies Cox — ^x^x and ujy 
ujyl \f\j are defined by 



and the position-momentum couphng is 



E 



(B6) 



(B7) 



To be consistent with the adiabatic limit, we can set 
A. - 1. 



The adiabatic elimination of the electronic excited levels 
leads to the AC stark shift 



^int 

where 



[Ai |5i> 



A2|52>(52|]cos2(^oa;), (C3) 



Ai = 







2 









4(a;sw - wi)' 



(C4) 



1 




2 




r _ 




— CJ2 



Since Ait)^, A27)^ <C \^y\ {fj = Kq^x < 1 is the Lamb- 
Dicke parameter), we can project i^int onto the ground 
state 11)^ oi (Jx^ and obtain 



i^int = ^int COS^{KQrx), 



(C5) 



where K^d = 47rn (n is an integer), ^lint is the effec- 
tive Rabi frequency. By assuming that the center of the 
microtrap is at the anti-node of the standing wave and 
expanding the cosine term in i^int, we obtain the phonon- 
phonon interaction 



H\j = Uintdx o^x^ 



(C6) 



where we have neglected the phonon number non- 
conserving terms. Notice that the second order in Vx just 
renormalizes the trap frequency, and can be absorbed in 
ujx' Here, the on-site interaction Ui^t = '^^intV^ can be 
tuned by the external laser beams. 
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